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ABSTRACT 

The Planck-ATCA Co-eval Observations (PACO) have provided flux density mea- 
surements of well defined samples of AT20G radio sources at frequencies below and 
overlapping with Planck frequency bands, almost simultaneously with Planck observa- 
tions. We have observed with the Australia Telescope Compact Array (ATCA) a total 
of 482 sources in the frequency range between 4.5 and 40 GHz in the period between 
July 2009 and August 2010. Several sources were observed more than once. In this pa- 
per we present the aims of the project, the selection criteria, and the observation and 
data reduction procedures. We also discuss the data in total intensity for a complete 
sample of 189 sources with 52oghz > 500 mJy, Galactic latitude |&| > 5°, and declina- 
tion 6 < —30°, and some statistical analysis of the spectral behaviour and variability 
of this sample, referred to as the "bright PACO sample" . Finally we discuss how these 
data could be used to transfer absolute calibrations to ground based telescopes using 
the CMB dipole calibrated flux densities measured by the Planck satellite, and we 
provide some test fluxes on bright calibrators. 
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1 INTRODUCTION 

The ESA's Planck satellit«E] is surveying the sky in nine fre- 
quency bands (30, 44, 70 GHz for the Low Frequency Instru- 
ment, LFI, and 100, 143, 217, 353, 545, 857 GHz for the High 
Frequency Instrument, HFI, with the beam FWHM ranging 
from 33 to 5 arcmin, Planck Collaboration 2011c). At the 
HFI frequencies it will provide the first all-sky surveys ever, 
while at LFI frequencies its higher sensitivity and resolution 
will allow a significant improvement over WMAP (Gold et 
al. 2011, Planck Collaboration 2006 and 2011c, Leach et 
al. 2008, Massardi et al. 2009). Planck thus offers a unique 
opportunity to carry out an unbiased investigation of the 
spectral properties of radio sources in a poorly explored fre- 
quency range, partially unaccessible from the ground. 

Planck's science yields can be greatly increased by si- 
multaneous (i.e. not affected by variability) ground based 
observations at lower frequencies as well as at frequencies 
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overlapping with Planck channels. Simultaneous observa- 
tions at overlapping frequencies are important to transfer 
the Planck flux density calibration, which is based on the 
CMB dipole, to the primary calibrators used by the ground- 
based radio telescopes. These observations can also be used 
to help with the validation of sources detected by Planck. 

At present, the flux density scales for ground-based tele- 
scopes are well determined for frequencies up to 10 GHz but 
there are known problems at the higher frequencies. The dif- 
flculty occurs because the radio interferometers, which can 
make the most precise flux density measurements at high 
frequencies, need point sources to calibrate and the strong 
point sources are mostly highly variable AGNs. This has 
led to the use of the planets. Mars and Uranus, as primary 
calibrators. These planets are resolved at all but the short- 
est interferometer baselines at the higher frequencies and 
their flux densities are estimated using relatively complex 
atmospheric models. In the case of Uranus there is a fur- 
ther complication because of the known orbital change in 
flux density as our view changes from equatorial to polar 
during the Uranian year. Simultaneous observations of very 
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bright sources give a unique opportunity to tie all the high 
frequency flux density scales together with the Planck's cal- 
ibration. It should be noted however that the source flux 
densities given in the Early Release Compact Source Cat- 
alogue (ERCSC, Planck collaboration 2011b) are averages 
over all observations, that may extend over several months. 
Therefore they may be used for some statistical purposes, 
as in Planck Collaboration (2011d) and in papers in prepa- 
rations from our group, while for calibration purposes we 
need to use the Planck time ordered data, averaged over 
shorter periods, comparable with the timescales of ground 
based observations. 

Source confusion is a serious issue for Planck, due to 
its rather large beam sizes. Sources near Planck's detection 
limit may be strongly affected by the Eddington's (1913) 
bias that leads to systematic flux density overestimates 
(Hogg & Turner 1998). Also, confusion fluctuations may pro- 
duce high intensity peaks that may be misinterpreted as real 
sources. Moreover, confusion can shift the positions of inten- 
sity peaks from the true source positions (Hogg 2001), com- 
plicating their identiflcation. Simultaneous ground based ob- 
servations, with much better resolution and signal-to-noise 
ratio will allow an accurate control of these effects. 

Source Spectral Energy Distributions (SEDs) over a fre- 
quency range as large as possible are crucial in determining 
their physical properties and in identifying the different com- 
ponents that may contribute to their emission. Given that 
observations in the full Planck frequency range will not be 
repeated at least in the foreseeable future, it is essential not 
to lose the opportunity for near simultaneous observations 
while Planck is flying. 

Observations simultaneous with Planck can be planned 
only for previously known sources. The Australia Telescope 
20 GHz (AT20G; Murphy et al. 2010; Massardi et al. 2010) 
survey provides the best sample for this purpose. This is 
the largest ground-based sample of the high radio frequency 
(> 10 GHz) sky: it is 93% complete above 100 mjy at 20 
GHz over the whole Southern sky with follow-up (within 
a few weeks time) at 4.8 and 8.6 GHz. For typical source 
spectra (Massardi et al. 2008), it is somewhat deeper than 
Planck, that is expected to reach detection limits ranging 
from ~ 400 mJy at 30 GHz to ~ 220 mJy at 100 GHz (Leach 
et al. 2008). 

The above considerations have motivated the Planck- 
ATCA Co-eval Observations (PAGO) project, described in 
§[2l The observations exploited the capabilities of the new 
Australia Telescope Compact Array Broadband Backend 
(CABB; Ferris & Wilson 2002, Koribalski et al. in prep.) 
system between 4.5 and 40 GHz. 



2 THE PACO PROJECT 
2.1 Sample selection 

The PACO sample is a combination of several sub-samples 
selected in order to maximize the scientiflc outcome in the 
allocated time. The ATCA allows us to make observations 
at frequencies lower than or close to those of the lowest fre- 
quency Planck channels at 30 and 44 GHz. We have ob- 
served 482 AT20G extragalactic sources at Galactic lati- 
tude |b| > 5° and outside a 5° radius circle centered at 




Figure 1. Polar projection of the Southern sky, showing the dis- 
tribution of the PACO samples (RA=Oh is at the top and the RA 
increases counterclockwise). The dotted lines indicate the Galac- 
tic plane and bound the area with Galactic latitude |6| < 5°. 
Sources in these regions were excluded from the PACO samples. 
The faint PACO sample (open pink diamonds) occupies the area 
with 3h < RA < 9h. The bright PACO sample (filled blue dia- 
mons) covers the whole area at 5 < —30° . The spectrally-selected 
sample (green squares) , the blazars and the AT calibrator sources 
(black asterisks) are spread over the whole southern sky. 



RA = 80.8939°, 5 = -69.7561° (the LMC region). Of 
these, 344 sources form 3 partially overlapping complete sub- 
samples, selected for different purposes. 

• The "faint PACO sample" is made of the 159 sources 
with S'20GHz > 200 mJy in the southern ecliptic pole re- 
gion (ecliptic latitude < —75°) and with 3h < RA < 9h, 
5 < —30° (Bonavera et al. 2011). Near the ecliptic poles 
Planck's scan circles intersect. Therefore the area is covered 
many times, and Planck's sensitivity is maximal. Neverthe- 
less, we expect that many sources in this sample will not 
be detected by Planck. On the other hand, Poisson fluctua- 
tions due to extragalactic sources below the detection limit 
(expected to be of 100-250 mJy in the range 100-353 GHz, 
Leach et al. 2008) are the main contaminant of Planck maps 
on angular scales smaller than 20'-30' at frequencies of up 
to ~ 100 GHz (Toffolatti et al. 1998, 1999; De Zotti et al. 
1999). It is therefore important to characterize as accurately 
as possible radio sources down to below the Planck detection 
limits at the time of Planck observations. 

• The "spectrally-selected PACO sample" comprises the 
69 sources with ^oghz > 200 mJy and spectra classified 
by Massardi et al. (2010) as inverted or upturning in the 
frequency range 5-20 GHz, selected over the whole southern 
sky. 
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Table 1 . The flux density values of the model for the primary cal- 
ibrator PKS B1934-638 at the center of our observing frequency 
bands, as it is in the Miriad version we used. 



Frequency 


Flux density 


[GHz] 


[mjy] 


5.5 


4965 


9.0 


2701 


18.0 


1101 


24.0 


725 


33.0 


441 


39.0 


335 



• The "bright PACO sample" comprising the 189 sources 
with S'20GHz > 500 mJy at 5 < —30°, is presented and dis- 
cussed in this paper. 

• An additional sample (partially overlapping with those 
defined above) is made of 63 bright southern blazars in- 
cluded in a long-term (up to 15 years) monitoring program 
at various frequencies with the Swedish-ESO Submillime- 
tre Telescope (SEST, Tornikoski et al. 1996) and observed 
also with the LABOCA instrument on the APEX. Based on 
a source counts analysis, this sample appears to be 91 per 
cent complete down to the cutoff of 1 Jy at 90 GHz. These 
sources are of special interest because some of them may be 
flaring during the Planck mission. If this happens, we get 
an unusually broad simultaneous frequency coverage of the 
flare. 

The samples described above include 76 ATCA calibrators 
with 20 GHz flux density S'2oghz > 200 mJy that showed 
more than 10% variability at this frequency in the last three 
years. Further 121 variable ATCA calibrators were added to 
the PACO sample, so that their total observed number is 
of 197. Our data will also be included in the AT calibrator 
database. 

Our observations allow us to characterize the source 
spectra in the range 4.5-40 GHz. Planck data will allow us 
to extend the spectral coverage up to 857 GHz, at least for 
the brightest flat-spectrum sources. 

The classification as extragalactic objects and an indi- 
cation of source extendedness up to 20 GHz were taken from 
the AT20G catalogue. AT20G positions were also used un- 
less more precise coordinates were available, as is the case, 
e.g., for the AT calibrators whose positions are based on 
VLBI observations. A map of our sample is shown in Fig. 
[T] As pointed out by Massardi et al. (2010) the sample may 
still contain some Galactic sources, but they should be very 
few because the relatively high resolution 20 GHz obser- 
vations preferentially select compact objects (< 30 arcsec) 
and Galactic sources with |6| > 1.5° are rarely that small, 
especially at the brighter fiux density levels. 

2.2 Observations 

The Planck scanning strategy has been publicly released 
soon after the satellite launch. The satellite scans the sky in 
circles passing close to the ecliptic poles. All Planck receivers 
cross the position of sources lying on the ecliptic equator 
within a few days. Sources close to the ecliptic poles remain 
within the satellite focal plane for up to 2.5 months. The 
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Figure 2. Model for PKS 1934-638 in the frequency range be- 
tween 4.5 and 40 GHz, as it is in the Miriad version we used. 



satellite covers the whole sky in about 6 months. Because the 
maps are the result of averaging all the scans through a po- 
sition, we can't get information from them on the very short 
term variability. Variability tests on the time ordered data 
are limited to the brightest objects because of lack of integra- 
tion on source. We consider observations to be 'co-eval' with 
the satellite if performed within 10 days from the satellite 
observations at any of its LFI frequencies. The Planck-On- 
the Flight forecaster (Massardi & Burigana 2010) has been 
developed and used to predict when the PACO sources are 
observed by the satellite, according to its pre-programmed 
pointing lists and the focal plane properties. 

The new ATCA CABB system allows 2x2 GHz simul- 
taneous bands in continuum. Applied to the 6 x 22 m ATCA 
antennas it gives a noise level of 0.5 mJy in 1 min on source 
in the 7 mm band. We choose to use the 7 mm receivers 
with bands centered at 33 and 39 GHz, to overlap the lower- 
frequency Planck channels, the 12 mm receiver with bands 
centered at 18 and 24 GHz, to include the AT20G selection 
frequency, and we extended the SEDs to lower frequencies 
using the 3-6 cm receiver with bands at 5.5 and 9 GHz. 

The AT20G catalogue is dominated by compact flat 
spectrum sources so the AT20G selection guarantees that 
most radio sources have small angular sizes so that the par- 
ticular configuration of the ATCA is not crucial for flux den- 
sity estimation. Hybrid arrays are preferable for equatorial 
sources, for which the beam shape of EW arrays is elongated. 
Because of the risk of confusion with the highly elongated 
synthesized beam for sources with declination 5 > —10°, 
they have not been observed at 3-6 cm. This is not a con- 
cern at shorter wavelengths where there is no confusion in 
the smaller primary beam. 

The scheduling process tries to guarantee the simultane- 
ity with the Planck observations following the satellite scan- 
ning strategy and to minimize the slewing time, maximize 
the time on source and complete the measurements at the 
selected frequencies (i.e. including considerations about the 
weather, the array conflguration in use and the minimization 
of the band switches). 

For each epoch and each frequency, a very bright com- 
pact object is used for bandpass calibration. Observations 
of the bandpass are usually repeated several times to ac- 
count for time instability of the band. PKS B1934-638 is 
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Table 2. The 63 good weather observational runs for the PACO project between July 2009 and August 2010. 
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4.5-25 GHz 


2010-08-14 


H168 


6.5 
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10 
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2010-01-27 


6A 


4 


4.5-25 GHz 


2010-08-21 


H168 


10 


4.5-40 GHz 


2010-02-06 


6A 


12 


4.5-25 GHz 











the primary calibrator for all the frequencies. The spectral 
model for the primary cahbrator that we have used in all 
the epochs is in Fig. [21 Tab. [1] lists the model flux density 
values at the center of our observing frequency bands. 

At 7 mm we have sometimes observed Uranus. Although 
the planet is always resolved at our longer baselines, we can 
use the shorter baselines to transfer the calibration scale. By 
comparing with the PKS 1934-638 flux densities, we noted 
that the model for Uranus requires a flux density correction 
by a factor of ~ 0.82 at 33 GHz and ~ 0.84 at 39 GHz. 
This is in agreement with what estimated on the basis of 
the on-going calibrator model update, as the absolute flux 
density scale at this frequency is still under investigation 
at both ATCA and VLA (J. Stevens, priv. comm., see Q. 
Hence we applied these correction factors to all the 7 mm 
flux densities calibrated with the planet. 

A leakage calibrator selected among the polarized ob- 
jects of the AT20G (a bright point-like unpolarized one 
should work as well) is observed at least 3 times, better if 
across the transit of the source, at each frequency for almost 
all the observing runs. Pointing calibration is performed reg- 
ularly for 12 and 7 mm observations. Each source at each 
frequency is observed for 1.5 min in a single pointing. 

Thanks to the compactness and high signal to noise 
of most of the PACO sources neither imaging nor phase 
calibration is necessary. Since sources are self-calibrated we 



get good images when observing with hybrid arrays. We got 
a suitable flux density estimation from the visibilities using 
the triple-correlation techniques, as will be explained in the 
next sections. 

The project got ~450 hours allocated over 65 epochs 
between July 2009 and August 2010 (i.e. during the time 
when Planck completed two all-sky surveys), of which 2 
epochs have been discarded for bad weather. Because of the 
high sensitivity in total intensity and the brightness of the 
sources in our sample we could observe on average 20 ob- 
jects per run at all the 3 pairs of frequencies. Bad weather 
seriously affected some of the Austral summer epochs mak- 
ing it impossible to get the 7 mm data in a few runs. Table[5] 
summarizes the configurations, and frequencies observed in 
each run. 



2.3 Data reduction 

Data reduction is performed by a C-shell coded pipeline 
which uses taskflfrom the MIRIAD software package (Sault 
et al. 1995) according to the scheme in Fig. [S] Chunks of 

^ A comprehensive description of each 

MIRIAD task can be found at the website 
jhttp: / / www.atnf .csiro. au / computing/ software / miriaxi / taskindex.html 
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Figure 3. Data reduction pipeline diagram. In capital letters are 
the Miriad tasks used. 



data affected by bad weather conditions or instrumental 
maffunction are removed, as well as data affected by shad- 
owing or known radio interference. Automatic procedures 
identify and flag among the 2048 channels, which constitute 
the 2.048 GHz bandwidth for each IF, those that exceed by 
5(T the mean value for each baseline for each source. Typi- 
cally less than 1% of the band is flagged. In order to properly 
deflne the detailed source spectral behaviour we have split 
each 2 GHz band into 4x 512 MHz sub-bands, and calibrated 
each sub-band independently. The steps for calibration are 
plotted in Fig. H) 



2.3.1 Extended sources 

We expect that the AT20G analysis of extended sources have 
already identified all the extended sources up to 20GHz. We 
have therefore done the analysis only for the 7 mm data, 
applying the automatic procedure to distinguish between 
point-like and extended sources used by Massardi et al. 
(2008) for the AT20G bright sample. 

The ratio between the actual and the theoretical phase 
closure rms calculated accounting for the noise in the beam 
gives an indication on whether the source is extended. For 
a point source the phase closure rms should be zero and is 
independent of instrumental effects. When the ratio of actual 
to theoretical phase closure rms is larger than 3 the source 
should be considered extended. We have used this ratio as 
an indication of extendedness for the sources that appear to 
be extended at 7 mm and were not identified as extended in 
the AT20G catalogue. 



2.3.2 Flux density esttmation 

To estimate the flux densities for point sources we have used 
the triple product (i.e. the geometrical average of the am- 




Figure 4. Scheme of the calibration performed with MIRIAD 
tasks. The indications in parenthesis refer to the case of Uranus 
as primary calibrator. The steps included in the dotted square 
have been applied only when a suitable leakage calibrator has 
been observed. 



plitude over a closure triangle, averaged over all the possi- 
ble triangles, given by the MIRIAD task CALRED). Our 
flux density estimation techniques are well suited for point 
sources. We adopted the amplitude of the shortest spac- 
ings for the array configuration used (see Tables [2] and |3]) , 
given by the MIRIAD task UVFLUX, as the best flux esti- 
mator. This approach still underestimates the flux densities 
for sources more extended than 15(100m/timin)(39GHz/!/) 
arcsec, where fe^in is the shortest baseline for any used con- 
figuration. Thus our fiux densities are increasingly underes- 
timated and, as a consequence, our samples are increasingly 
incomplete for more and more extended sources. Table[3]lists 
the shortest baselines (&min), and the primary and synthe- 
sized beam sizes for the central frequencies of our observing 
runs for all the array configurations used. 

The noise term of the fiux density errors for point 
sources are calculated as the rms of the triple product am- 
plitudes over all the possible triangles of antennas. For ex- 
tended sources the rms of the triple product has been mul- 
tiplied by ^/r^ where ni, is the number of baselines of the 
used configuration. An analysis of the gain term of the fiux 
density error including systematics and calibration uncer- 
tainties is presented in ij |2.4.1l 



2.4 Quality assessment of the results 

By running the data reduction pipeline previously described 
we have generated a preliminary catalogue for each observ- 
ing epoch. This catalogue was then processed through a 
quality assessment procedure. Depending on the results of 
the assessment, we either proceeded to create the final cat- 
alogue for the considered epoch or we carried out further 
investigations in order to improve the results. 
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Table 3. The array configurations and their resolution properties. 



Array 


Shortest 


frequency 


Primary beam 


Beam IFl 


Beam IF2 


^max 


Configuration 


BaseUne [m] 


IFl - IF2[GHz] 


FWHM [ 


arcmin] 


[arcsec] 


[arcsec] 


[arcmin] 


H75c 




5.5 - 9 


8.7 - 


5.3 


1 1 ^ ^1 V 1 1 ^ ^1 

J- J- U . U J- jV J.J_0.tJJ- 


tu.oij X lU.oy 


6.0 - 


3.7 






18 - 24 


2.6 - 


2.0 




9fi 47 Y 9fi 47 


1.8 - 


1.4 






33 - 39 


1.5 - 


1.2 


1 Q 9^ V 1 Q 9^ 


1 fi 9Q Y 1 9Q 


1.0 - 


0.9 


H168d 


U J- I -L 


5.5 - 9 


8.7 - 


5.3 


^fl 9*^ V ^0 9"^ 


•^0 fiQ X "^0 fiQ 


3.1 - 


1.9 






18 - 24 


2.6 - 


2.0 


l-O.OO J'^ ±0.00 


11 ^1 IT 1 1 ^1 


0.9 - 


0.7 






oo - oy 


1.5 - 


1.2 


O.O I X o.o ( 


7 n^i Y 7 OR 


0.5 - 


0.4 


H214c 


82 (2-4) 


5.5 - 9 


8.7 - 


5.3 




23 78 X 23 78 


2.3 - 


1.4 






18 - 24 


2.6 - 


2.0 


11 J^Q -if 1 1 SQ 


Q9 Y 8 Q9 
o. y ^ A- o . 


0.7 - 


0.5 






33 - 39 


1.5 - 


1.2 


6.48 X 6.48 


5.49 X 5.49 


0.4 - 
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The quality assessment step is needed because each ob- 
servation can be affected by different problems (bad weather 
conditions, some antenna correlations to be rejected, resid- 
ual spikes, etc.) which could not be addressed by the auto- 
matic pipeline. It is designed to identify the major ffaws the 
catalogue of each observing epoch may have and to perform 
an a-posteriori flagging of the bad data. 

The quality check assumes that the true SEDs do not 
contain discontinuities; those observed are interpreted as due 
to residual spikes in the data or to calibration errors. The 
smoothness of the spectra is checked by comparing the data 
points with a fitting polynomial function with degree equal 
to the number of 2 GHz bands observed. We stress that, at 
this stage, the fit is done for flagging purposes only. 

For each SED of the considered epoch we compute the 
fractional rms divergence, cr, of the data points from the 
fitting line. Next we exclude those points (if any) that are 
more than 4a off and recompute the fractional divergence 
over all the other points. The procedure is iterated until the 
value of a converges. To get rid of particularly bad epochs, 
if the final a is higher than 0.02, which corresponds to ~ 2x 
the median rms fractional divergence for all the observing 
epochs, we do one more iteration with a — 0.02. Whenever 
two data points within the same 2 GHz band are discarded, 
it is likely that some problems are affecting the whole band, 
and we flag it for the considered source. 

When the fraction of good data for a given observing 
epoch exceeds 90% (and this happens in most cases) we di- 
rectly create the final catalogue of the unflagged data. All 
the other cases are inspected in order to identify the prob- 
lems and new runs of the data reduction pipeline with im- 
proved specifications are performed aiming at reducing the 
data loss to no more than 10-15%. The recovery succeeded 
in most cases. The exceptions are due to problems affecting 
observations of the flux density calibrator for some bands. 
In a few cases the calibrator couldn't be recovered to a satis- 



factory level and the whole target sample was rejected. After 
the quality assessment process, we are left with ~ 85% of 
our data. 



2.4-1 Flux density error estimation 

The gain term of the noise for each epoch was estimated as 
the rms fractional divergence between the data and the flt, 
re-calculated as discussed in the previous section after all 
the flagged data have been removed. For the days with less 
than 10 good flts available we have used ag — 0.012, equal to 
the median of the rms fractional divergence calculated over 
all the epochs for point sources only. The gain term, as we 
calculated it, includes also flux density errors due to pointing 
errors. Hence, the final error bars for each point-source in 
each epoch have been estimated as square root of the square 
of the gain term multiplied by the source fiux density plus 
the square of noise term (estimated as discussed in I2.3.2p . 
Since our flux density estimation techniques are not well 
suited for extended sources, we have increased their gain 
term used in the error estimate to a minimum of ~ 5 per 
cent of their flux density. This value has been obtained by 
scaling the rms fractional divergence calculated over all the 
epochs by where rb = 15/1 is the ratio between the 
number of baselines used in case of point sources and in 
case of extended sources. This term does not account for 
any missing flux density of extended objects. 

The gain term does not include any scaling error due to 
deficiencies in the absolute calibration errors which should 
affect equally all the frequencies calibrated with the same 
calibrator source (generally PKS 1934-638) . The comparison 
with Planck CMB-dipole calibrated flux densities will allow 
us to estimate this bias and correct it. 
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Figure 5. Spectra of sources in the PACO bright sample (crosses). Different colours correspond to different observing epochs. Diamonds 
corresponds to epochs simultaneous to Planck observations. To avoid overcrowding we have not plotted the error bars that are in all 
cases smaller than the symbol size. Also shown, for comparison, are the AT20G data (triangles) and, when available, WMAP data taken 
from the NEWPS_5yr catalogue (squares; Massardi et al. 2009). 



M. Massardi et al. 



CO 1 .0 
0.1 



M20GJ05591 1 -452940 
^* 


M20GJ060031 -393701 

"^m — 


M20GJ050425-422530 
*^ % 


AT20GJ062307-643620 

*, 

>^ * X W 


M20GJ062706-352916 


M20GJ0631 1 1 -415427 


M20GJ063546-751616 

mm 


M20GJ06481 4-30441 9 

'^M * 


M20GJ0701 34-463436 

HH ;|| ifl ^ 


AT20GJ074 145-470926 


CT20GJ074320-561933 


M20GJ074331 -672625 


M20GJ081 108-492943 


\T20GJ082525-501039 


M20GJ083529-59531 1 
m **** 


AT20GJ084502-545808 
^ttf ^ Mm 

>jgK ^ ^ 


M20GJ084945-354101 


M20GJ090453-573504 


M20GJ092246-395935- 

mil JH 


M20GJ 1001 59-443800 


M20GJ1 00731 -333306 
* * * 


AT20GJ101450-450841 


M20GJ101828-312353 

mm mm 


M20GJ102343-664647 



10 10 10 10 10 10 

freq[GHz] freq[GHz] freq[GHz] froq[GHz] freq[GHz] freq[GHz] 



10.0 

ui 1 .0 

0.1 
10.0 

CO 1 ,0 

0.1 
10.0 

1 .0 

0.1 
10.0 

■Tl 1 .0 

0.1 



M20GJ103653-374415 

mw 


\T20GJ 104 144-473959 


M20GJ1 051 09-534446 


AT2GGJ 105843-800353 

to* iiit 


M20GJ1 10204-440422 


AT20GJ1 10352-535700 

mW * ^ ««■ 


M20GJ1 10708-444907 

E 


\T20GJ1 10712-682051 


M20GJ1 1 1827-463415 

^ Ml ^ 


AT20GJ1 13601-682705 


M20GJ1 14531-483610 

* 

* 


AT20GJ1 14553-695404 


M20GJ1 14701 -381 210 


AT20GJ1 14733-675341 


M20GJ1 15253-834410 


AT20GJ1 15421 -350529 

% 


M20GJ121806-460030 


AT20GJ122454-S31310 


M20GJ122726-443639 
5«<« X*X)C 


M20GJ 124828-455947 


M20GJ125459-713818 


AT20GJ1 25759-31 551 5 

^^^^ 


M20GJ1 30349-554031 


AT20GJ131504-533436 



10 10 '0 10 10 10 

freq[GHz] freq[GHz] freq[GHz] freq[GHz] freq[GHz] freq[GHz] 



Figure 5. Continued. 
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Figure 5. Continued. 
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Figure 5. Continued. 



3 THE BRIGHT SOURCES PACO 
CATALOGUE 

Of the 189 sources comprised in the PACO bright sample, 13 
are flagged as extended in the AT20G catalogue and 2 more 
(J062706-352916 and J130527-492804) were found to be ex- 
tended by our 7mm data analysis. Two of the 13 AT20G ex- 
tended sources (J051926-454554 and J051949-454643) cor- 
respond to the core and to the western lobe of Pictor A. 
We have performed a dedicated observation in total inten- 
sity and polarization at 18, 24, 33 and 39 GHz of this well- 



known extended radio source by mosaicking the region in- 
cluding the core and the western lobe, which is the brightest 
and the more highly polarized (Burke-Spolaor et al. 2009). 
These observations will be presented in a separate paper and 
the two sources have been left aside. 

After the data reduction, 7 extended sources (J040848- 
750720, J084127-754028, J130527-492804, J161505-605427, 
J193557-462043, J195817-550923,J235904-605503) turned 
out to have unreliable flux density estimates at all the fre- 
quencies and for all the days of observations. This is mostly 
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Figure 6. Images of the source AT20GJ053850-440508, observed 
on July 26th, 2009 with the hybrid configuration H75, as an ex- 
ample of the quality of our data and a representation of the beam. 
Each panel enclosing an image has a ^ 6' X 6' size. Rows show 
the images obtained in the four 512-MHz sub-bands associated to 
each of the 6 IFs. 



due to having used too extended array configurations, so 
that the sensitivity on the observed scale was poor. They 
have therefore not been included in the present catalogue. 

For all the other extended objects (J052257-362730, 
J062706-352916, J074331-672625, J114531-483610, J132527- 
430104, J215706-694123) the estimated flux densities are 
lower limits to the integrated flux densities. Since these 
sources passed all the quality checks they have been included 
in the catalogue, with flux densities listed as lower limits, but 
have been excluded from further analysis. 

Because of the large size of the table, the catalogue of 
the 180 sources for which we have a flux density estimation 
or a lower limit is available only online. Columns are as 
follows: 

1 AT20G name; 

2-3 J2000 equatorial coordinates; 

4 epoch of observation as [yymmdd] ; 

5 flag 's' for "simultaneity" (i.e. within 10 days from the 
Planck observations). The epochs before February 2010 have 
to be compared with the flrst Planck survey, whereas the 
following epochs should be compared with the second Planck 
survey. 

6 flag 'e' for extended sources; 

7-31 flux densities in mjy at the 24 frequencies between 
4732 and 39768 MHz in which our 6x2-GHz bands have 
been divided; 

32-55 flux density errors; 
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Figure 7. Colour-colour plot comparing the spectral indices in 
the ranges 5-10 and 30-40 GHz. The dotted lines shows the bound- 
aries adopted for our spectral classification. 



56-60 flt parameters, ^o, fo, o,, b (see 8^ 13. Ijl . 

Figure [5] shows the spectra for the sources in the PACO 
bright sample. Different colours correspond to different ob- 
serving epochs. Figure |6] shows, as an example, images at 
each frequency for the point source AT20G J053850-440508. 
For the epochs observed with hybrid arrays, as in this case, 
good quality images can be obtained even if we observe only 
one cut per source. This is possible thanks to the bright- 
ness of the sources in our sample. It is a clear indication of 
the quality of our observations and of the beam size at the 
different frequencies. 



3.1 Spectral behaviour 

Even though a power-law is a poor description for the ma- 
jority of the spectra, conventional spectral indices are still 
useful for practical purposes. For example, they provide the 
easiest way to extrapolate the observed counts or model pre- 
dictions from one frequency to another. 

As our observations cover a wide frequency range (4.5- 
40 GHz), over which a single power-law is not enough to 
describe the spectral behaviour of the sources, we have stud- 
ied the spectra of our 174 point-like sources by fitting the 
observed data with a double power law: 



(1) 



where i/ is the frequency, 5* the fiux density in Jy and So, 
fo, a and b are free parameters. We considered only the 
sources for which we had at least 4 data points for each of 
the 2x2 GHz bands considered. The fit has been performed 
in logarithmic units by minimizing the with a nonlinear 
optimization technique based on an implementation of the 
Generalised Reduced Gradient optimisation method (GRG, 
Windward Technologies, 1997). The allowed ranges for the 
parameters are: 4.5 GHz ^ i/o 40 GHz, a —3, fe ^ 3. We 
have considered as bad fits those with x'^~{X^) > 3(7^ where 



ix ) = 1-3 and 



0.54 are mean and standard deviation 
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Figure 8. Distribution peak frequencies in the observer's frame 
(hatched histogram) and in the source frame (shaded) for sources 
with measured redshift. 



Figure 9. Distribution of the differences between low- and high- 
frequency spectral indices (al" — ctglj). 



of the Gaussian fit of the reduced distribution of the fits 
over the analyzed sample. 

For 4 of the 174 sources we do not have enough data 
to obtain meaningful values of the parameters, while in 5 
cases eq. ([ij does not yield an acceptable fit (x^ above the 
threshold) . For the remaining 165 sources we have classified 
the spectral shape according to the values of the low- and 
high-frequency spectral indices al^' and q:|o ('S' oc u") de- 
fined in the ranges 5-10 GHz and 30-40 GHz and calculated 
using the fitting formula. The distribution of these spectral 
indices is shown in Fig. [T] 

When — alol < 0.35 the source spectrum is classi- 
fied as 'single power law'. We have 24 such sources (~ 14%), 
further subdivided in: 

• 6 (3.6%) 'steep' with aj", a||] < -0.3, 

• 17 (10.3%) 'fiat' with \al°\, =^ 0.3, 

• 1 (0.6%) 'inverted' with al° , > 0.3. 

The other 141 have been classified as follows (see Fig. (Tj: 

• 'peaked' if al° > 0.3, < --0.3: 24 (14.5%) sources; 

• 'down-turning' if q|o ^ min(Q5'' — 0.35,-05") and 
< 0.3: 109 sources (66.0%); 

• 'self-absorbed' if 05° - 0.35 ^ af^ ^ max(-0.3, -al,°): 
8 sources (4.8%). 

Remarkably, we don't find any 'upturning' source (05" < 

030)- 

The peak frequency of 'peaked' sources is — 
VQ{—h/a)^^^''^~°'^\ The distribution of Vp in the observer's 
frame is shown by the solid histogram in Fig. [S] The mean 
value is 16.4 GHz with a standard deviation of 3.5 GHz. The 
median is 15.5 GHz with inter-quartile range 14.0-19.7 GHz. 
The figure also shows the distribution of peak frequencies in 
the source frame for the 9 sources for which redshifts are 
available (Massardi et al. 2010, Mahony et al. in prep.). For 
these the mean is 32.1 GHz with dispersion of 8.8 GHz; the 
median is 32.2 GHz with inter-quartile range is 27.7-38.9 
GHz. 

The distribution of the differences between low- and 
high-frequency spectral indices {01^ —a^) is shown in Fig. (9] 
There is a clear steepening of the spectra at high frequency. 
The mean difference is 0.75 with a standard deviation of 
0.37. The median difference is 0.73 with inter-quartile range 
0.52-0.97. 
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Figure 10. Distributions of spectral indices between 5 and 20 
GHz (solid line) and between 20 and 40 GHz (dashed line). 



Figure [TOI shows the distributions of spectral indices be- 
tween 5 and 20 GHz and between 20 and 40 GHz (we refer 
to 20 GHz since it is the selection frequency of the PAGO 
sample). The mean value of af is -0.07, with a standard 
deviation of 0.32, while the mean a|o is -0.55, with a stan- 
dard deviation of 0.34. It should be noted that the 20 GHz 
selection favours sources that are brighter at this frequency. 
Therefore a positive value of 01^ — a|o is, at least partly, 
due to a selection effect. 

Finally, Table |4] summarizes the median spectral indices 
calculated from the fitting double power-law function among 
couples of frequencies in the range covered by our observa- 
tions. Since spectra are complex, spectral indices were com- 
puted only for relatively narrow frequency intervals. 

Table 4. Matrix of median spectral indices calculated from the 
fitting double power-law function. 



Frequency [GHz] 10 15 20 30 40 



5 0.04 -0.02 -0.07 

10 - -0.12 -0.19 

15 - - -0.25 -0.38 

20 _ _ _ _o.45 -0.55 

30 _ _ _ _ _o.70 
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Table 5. Median variability index of the bright sample at differ- 
ent frequencies and time lags. 
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Figure 11. Top: median variability index of the sample as a func- 
tion of frequency; error bars are 1.2ba/\/ N — 1. Bottom: number 
of sources. 

3.2 Variability 

Most of the sources in the sample have been observed more 
than once to monitor the flux density variations between 
the first two Planck surveys, carried out six months apart. 
Thus we typically re-observed each source after 6 months. 
More frequent monitoring was carried out only if spare time 
was available. Therefore observations with time lags shorter 
or longer than 6 months were done only for limited sub- 
samples. 

Following Sadler et al. (2006) we define the variability 
index as: 

_ ioo./ E[^»-('g)]^-E'^? 

^'"^^ " (S) V n 

where Si is the fiux density of a given source measured at 
the i-th epoch, en is the associated error, n is the number 
of measurements, and {S) is the mean fiux density of the 
source, computed using all the available observations of the 
source (not only those used for computing the variability 
index) . 

To estimate the variability over a certain time interval 
we have selected for each source the best pair of observations 
(i.e., n = 2) spaced by the that time interval within ±20%. In 
Fig [11] we show the median variability index for the sample 
of sources as a function of frequency for time intervals of 3, 
6, and 9 months (upper panel) and the number of sources 
having pairs of observations for those time intervals (lower 
panel) . 

We find, as expected (Impey & Neugebauer 1988; Cia- 
ramella et al. 2004), a trend toward an increase of the vari- 



ability amplitude with frequency (see Tableland Fig. Illf) . 
There is also a marginal indication of higher variability for 
the longer time lag, consistent with earlier results (e.g. Cia- 
ramella et al. 2004). On the other hand, around 20 GHz the 
median variability index for a 9 months lag is ~ 9 ± 1.7%, 
slightly larger than that found by Sadler et al. (2006) over 
a one year timescale (6.9%) for a somewhat fainter sample 
(flux density limit of 100 mjy). The distribution of Kms, 
however, is similar to that found by Sadler et al.; for exam- 
ple, 37 ± 9% of our sources show a variability index larger 
than 10%, to be compared with 42 ± 7% for the Sadler et 
al. sample. 



4 ABSOLUTE CALIBRATION FOR GROUND 
BASED TELESCOPES 

The absolute calibration of the Planck detectors up to 353 
GHz is derived from the annual modulation of the CMB 
dipole by the satellite orbit around the Sun. An absolute 
measurement of the dipole is obtained by differentiating 
along a spin period. The CMB dipole direction is expected 
to be recovered to better than ~ 2 arcmin, and its ampli- 
tude to better than 0.5% (The Planck collaboration 2006). 
Planck aims at obtaining an absolute photometric calibra- 
tion better than 1% for all frequency channels up to 353 
GHz. As discussed in Lopez-Caniego et al. (2007) the main 
limitation to the accuracy of flux density estimates of sources 
extracted from low-resolution maps (appropriate for CMB 
studies), such as these produced by WMAP and Planck, is 
confusion due to faint sources within the beam. Leach et al. 
(2008) estimated typical rms errors ~ 100 mJy for the 30 
and 44 GHz channels at high Galactic latitude. This corre- 
sponds to a few percent of the fiux density of the brightest 
extragalactic sources. 

Thus ground based observations, simultaneous with 
Planck's and at approximately the same frequencies, of very 
bright sources can be used to assess the absolute flux den- 
sities of calibrators for ground based telescopes. The PACO 
project has enough observations to determine the relative 
flux density scale between Planck and the ATCA. 

Table[6]lists the 7 mm PACO flux densities of 6 compact 
sources with 20 GHz fiux density S'at20G > 5 Jy and high 
quality observations simultaneous to Planck's. With the ex- 
ception of AT20GJ053850-440508, these sources lie outside 
the region covered by the bright sample. All of them are 
well-known highly variable blazars. Two of them are equa- 
torial sources and are currently used in some ad-hoc ob- 
servational runs for cross-calibration between southern and 
northern hemisphere telescopes, in order to assure that the 
Planck absolute calibration could be shared with all the ma- 
jor ground-based facilities. 

As noted before, the ERCSC (Planck collaboration 
2011b) fluxes are unsuitable for calibrating ground-based 
telescopes. A proper comparison with Planck fiux density 
scales will be possible, and will be done, when Planck time 
ordered data will become available. A preliminary compar- 
ison between ground based and Planck measurements has 
been presented by Planck collaboration (2011b, 2011e). 
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Table 6. The 7 mm flux densities, measured simultaneously with Planck observations, of point-like calibrators in the PACO sample with 



»S'at20G > 5 Jy. The figure immediately below each flux density measurement 


is the corresponding error. 








AT20G name 


RA 

[hr] 


5 

[deg] 


Observing 
Date 


2GHz 
[mJy] 


^32 7GHz 

ImJy] 


^33 2GHz 

ImJy] 


^33 7GHz 

ImJy] 


Sgg 2GHz 
[mJy] 


^^38 7GIlz 
[mJy] 


Sgg 2GHz 

[mJy] 


'^39 7GHz 
[mJy] 










CTg32 2GHz 
[mJy] 


■^332 7GHz 
[mjy] 


c S33 2GHz 
[mjy] 


f 7GHz 
[mJy] 


'^SSS 2GHz 
[mJy] 


f Q Qta "71"' "Hr, 

[mJy] 


ooy . zijtiz 

[mJy] 


CQ ^lO "TC Hr, 

[mJy] 


J053850-440508 


5.6473226 


-44.0858154 


2010-02-28 


9506 


9483 


9385 


9356 


8962 


8921 


8864 


8893 










lis 


118 


117 


117 


112 


111 


110 


111 


,1053850-440508 


5.6473226 


-44.0858154 


2010-03-15 


14278 


14210 


14135 


14137 


14848 


14620 


14639 


14814 




















192 


192 


194 


J053850-440508 


5.6473226 


-44.0858154 


2009-09-15 


6344 


6352 


6339 


6339 


6192 


6156 


6177 


6141 










92 


92 


92 


92 


89 


89 


89 


89 


J125611-054721 


12.9364349 


-5.7893119 


2010-07-16 


14740 


14665 


14669 


14659 


14041 


13956 


14064 


14089 










126 


126 


126 


126 


120 








,1133739-125724 


13.6277171 


-12.9568596 


2010-07-16 


4228 


4197 


4188 


4181 


3974 


3955 


3922 


3876 










36 


36 


36 


35 


34 


34 


33 


33 


J183339-210341 


18.5610881 


-21.0611248 


2009-09-26 


3256 


3210 


3174 


3132 


2750 


2714 


2674 


2639 










16 


16 


16 


16 


14 


13 


13 


13 


J192451-291430 


19.4141825 


-29.2417011 


2010-04-01 


17111 


17807 


17451 


17748 






















196 


204 


200 


203 














,7192451-291430 


19.4141825 


-29.2417011 


2010-04-10 


14571 


14696 


14686 


14357 


14007 


13782 


13723 


13891 










213 


215 


215 


210 


205 


202 


201 


203 


J192451-291430 


19.4141825 


-29.2417011 


2009-09-26 


14028 


13876 


13789 


13663 


12490 


12393 


12263 


12127 










71 


71 


70 


69 


63 


63 


62 


62 


J192451-291430 


19.4141825 


-29.2417011 


2009-10-07 


15333 


15248 


15083 


14820 


14380 


14225 


14341 


14617 










60 


60 


59 


58 


56 


56 


56 


57 


,1192451-291430 


19.4141825 


-29.2417011 


2009-10-14 


13728 


13641 


13475 


13267 


12282 


12179 


12076 


12078 










122 


122 


120 


118 


110 


109 


108 


108 


J222547-045701 


22.4297933 


-4.9503870 


2009-11-06 


6177 


6098 


5929 


5905 


5557 


5337 


5009 


4978 










105 


103 


101 


100 


94 


90 


85 


84 


J222547-045701 


22.4297933 


-4.9503870 


2009-11-19 


5589 


5524 


5469 


5417 


4910 


4840 


4794 


4743 










110 


109 


108 


107 


97 


95 


94 


93 



5 CONCLUSIONS 

The PACO project consists in observations with ATCA of 
several complete samples in three pairs of 2 GHz bands (5.5 
and 9 GHz, 18 and 24 GHz, 33 and 39 GHz). The observa- 
tions were carried out in several epochs between July 2009 
and August 2010. At least one observation was made within 
10 days from the Planck satellite observations in any of the 
LFI channels (30, 44, 70 GHz). 

We have presented and analyzed the data on sources 
in the "bright sample" (S20GHZ > 500 mJy). Although the 
present analysis is preliminary, in wait for the much broader 
spectral coverage that will be possible combining our mea- 
surements with Planck data, we note some interesting re- 
sults: 

• Only a minor fraction (~ 14%) of source spectra have 
a shape close to a single power-law over the 4.5 to 40 GHz 
range; they are mostly "flat" -spectrum. 

• Few sources have a positive (i.e. rising, S oc u") high- 
frequency (30 to 40 GHz) spectral index, and none has an 
upturning spectrum (qsg'eS^ < o^mghz)- Most sources show 
a spectral steepening above 10-20 GHz: the mean value of 
the spectral index is ~ —0.07 (with a standard deviation 
of 0.35) between 5 and 20 GHz and decreases to ~ —0.57 
(with a standard deviation of 0.32) between 20 and 40 GHz. 
Break frequencies are almost uniformly distributed over the 
full frequency range covered by our observations. We do not 
find any correlation between the variation of the spectral 
index and the break frequency. We note that part of this 
steepening may be due to a selection effect since the primary 
20 GHz selection of the sample emphasizes sources which are 
brighter at this frequency. 

• A significant fraction (about 14%) of the sources show 
a peak within our frequency range. The mean value of the 
peak frequency is 16.3 GHz with a standard deviation of 3.5 
GHz. The mean rest frame peak frequency for the 9 such 
sources for which a redshift is available is 32.1 GHz with 



dispersion of 8.8 GHz. In general, however, the peak is not 
very prominent. 

• A handful of sources (8, i.e. 4.7%) show a rising low- 
frequency spectrum, followed by a flattening, perhaps in- 
dicative of self-absorption up to ~ 10 GHz. 

• There is a trend toward an increase of the variabil- 
ity amplitude with firequency and a marginal indication of 
higher variability for the longer time lag. 

We also present accurate ATCA flux density measure- 
ments, simultaneous with Planck, for 6 bright (S'at20G > 
5Jy) point-like sources that can be used to exploit the ex- 
cellent absolute calibration achieved by Planck to re-assess 
the high frequency calibration of ground-based telescopes 
and to improve models for planets and other sources used 
as calibrators at mm wavelengths, once the satellite time 
ordered data are available. 

On the other hand, our data will help in the validation 
of sources detected by Planck and in quantifying the effect 
of confusion on Planck flux density and position estimates. 
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